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ABSTRACT 

We evaluate the infinite volume, continuum limit of glueball masses 
in the valence (quenched) approximation to lattice QCD. For the lightest 
states with J^^ of 0++ and 2++, we obtain mo = 1340 ± 160 MeV and 
m2 = 1900 ± 320 MeV. 
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1 INTRODUCTION 



Although dozens of observed hadrons have been convincingly interpreted as 
bound states of combinations of quarks and antiquarks, physical particles composed 
primarily of chromoelectric field, glueballs, have so far not been identified unambigu- 
ously in experiment. One possible explanation is that glueballs simply may not exist. 
All excitations of the chromoelectric field might occur only mixed with a large num- 
ber of quark- ant iquark states. Each mixed state might have only a small probability 
amplitude of excited chromoelectric field and, as a result, show no easily identifiable 
experimental signature of the presence of a field excitation. Alternatively, some of 
the particles already observed in experiment may actually be glueballs and may not 
have been identified because QCD's predictions for glueball masses and decays have 
not yet been determined accurately enough. 

In the present article, we evaluate the infinite volume, continuum limit of the 
predictions of lattice QCD, in the valence approximation, for the masses of the lightest 
glueballs with J^^ of O^"*" and 2"*"+. The 0"'"+ is probably the lightest member of the 
glueball spectrum and the 2"*"+ the next lightest [|T], Q . The valence approximation may 
be viewed as replacing the momentum dependent color dielectric constant arising from 
quark- ant iquark vacuum polarization with its zero- momentum limit @]. For glueball 
masses, this approximation amounts to a reinterpretation of the predictions of pure 
gauge theory. Pure gauge QCD certainly predicts the existence of unmixed glueballs. 

The valence approximation could be expected to be fairly reliable for low-lying 
hadron masses, which are determined largely by the low momentum behavior of the 
chromoelectric field. This expectation is supported by a recent valence approximation 
calculation of the masses of eight low-lying baryons and quark-antiquark mesons. 
The predicted masses are all within 6% of experiment and all differences between pre- 
diction and experiment are consistent with the calculation's statistical uncertainties. 
The prediction we obtain here for the 0"*"^ glueball mass, 1340 ± 160 MeV, lies below 
the largest mass, = 1672 MeV, for which the valence approximation has now, to 
some degree, been established as reliable. Allowing a disagreement of up to 80 MeV, 
6% of 1340 MeV, between the valence approximation and full QCD, the predicted 
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mass is consistent with the masses of either of two well-estabhshed states, 
/°(1400) and /°(1590). The mass prediction for the 2++ glueball, 1900 ±320 MeV, is 
close enough to the range for which the valence approximation has been tested that 
it appears reasonable to expect a disagreement of at most 6% with full QCD in this 
case also, giving an additional uncertainty of 120 MeV. The overall uncertainty in 
the 2+"*" mass is large enough that this state could be any one of a half dozen or so 
observed resonances. 

Infinite volume, continuum limit glueball masses obtained from the valence 
approximation to lattice QCD have now also been reported in Ref . P| , which appeared 
while the present work was in progress. The predicted 0++ and 2++ masses are 
1550 ±50 MeV and 2270 ±100 MeV, respectively. The prediction for the O'*"^ mass is 
about 1.3 standard deviations above our result. It seems to us that the central values 
of both predictions may be somewhat larger and the error bars of both predictions 
somewhat smaller than the data reported in Ref. actually suggests. We will return 
to this subject below. 

The calculations described here were done on the GFll parallel computer 
at IBM Research [|| and took about six months to complete. GFll was used in 
configurations ranging from 375 to 480 processors, with sustained speeds ranging 
from 5 Gflops to 7 Gflops. 

2 DEFINITIONS AND ALGORITHMS 

To define propagators from which light glueball masses can be extracted effi- 
ciently, it is useful to form smeared gauge fields with reduced coupling to the under- 
lying gauge field's high momentum components 0, § . In place of the gauge invariant 
smearing of Refs. 0, H, however, we use a smearing method based on Coulomb 
gauge. Each gauge configuration on which glueball propagators are to be constructed 
is transformed to lattice Coulomb gauge, defined to maximize at each site the target 
function X^xi -ReTr[f/j(x)], where the sum is over all lattice sites and space direction 
gauge links. A gauge transformation which produces a local maximum of this sum is 
found by a method qualitatively similar to the Cabbibo-Marinari-Okawa Monte Carlo 
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algorithm. The lattice is swept repeatedly, and at each site we maximize the target 
function first by a gauge transformation in the SU (2) subgroup of SU (3) acting only 
on gauge index values 1 and 2, then by a gauge transformation in the SU (2) subgroup 
acting only on index values 2 and 3, then by a gauge transformation in the subgroup 
acting only on index values 1 and 3. Maximizing the target function over SU{2) 
subgroups is easier to program than a direct maximization over all of SU{3). On the 
other hand, it is not clear that maximizing each site over SU{3) would significantly 
accelerate the full transformation to Coulomb gauge. A local maximum is reached 
when at each site the quantity R{x) vanishes where 

Q{x) = Y.iU^{^)-UK^-^)]: 

i 

R{x) = Q{x)-Q\x)-\lmTr[Q{x)]. (2.1) 

The vector i is a unit lattice vector in the positive i direction. We stop the iteration 
process when the sum over the lattice of the quantity Tr[R\x)R{x)] becames smaller 
than a convergence parameter r. 

Coulomb gauge smeared fink variables U?{x) for s > 1 are then defined by the 
average 

^^(^) = 77TT^ E U,{x + pj + qk), (2.2) 

V* "I" ^) 0<p,q<s 

where j and k are the two space directions orthogonal to i. The product of s sequential 
Utiy) gives Vi{x), 

Vax)^Ut{x)...Ut[x + {s-l% (2.3) 

From the V^{y) we define smeared loops L^j{x) 

Ll.{x) = Tr[Vl{x + si)V^'{x + 2si + sj)y/(x + si + 2s]) V/(x + sj% (2.4) 

This definition fits together as smoothly as possible the products of smeared links 
forming each side of Lfj{x). The sum of L^j{x) over all a; in a hyperplane with fixed 
time component t defines the zero- momentum loop variable Lhit). 
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For smearing size 0, the loop variable L^j{x) is defined to be an unsmeared 
plaquette 

L'',.(x) =Tr[U,{x)Uj{x + t)U,{x +jyU,{x)^], (2.5) 

and L^j{x) is then defined from L^ij^x). 

A field for the 0+"'" glueball, transforming according to the Ai representation 
of the cubic symmetry group [0], is given by 

A\t) = Y.Re[UM- (2-6) 

Fields for the 2"*""^ glueball, transforming according to the E representation of the 
cubic symmetry group, are given by 

El{t) = Re[L{,{t) - LUt)l 

Em = Re[L{m + Lim-2L\,{t)]l^ (2-7) 
From these glueball fields, we define glueball propagators as 

Cf{t) = < A'{t)A''{0) >-<A'{t)>< A'\0) >, 

Cr'it) = Y.[<E-it)EfiO)>-<E^it)><Ef{0)>], (2.8) 

i 

Although < E-{t) > vanishes for a lattice which is invarient under the cubic symmetry 
group acting on the space directions, one of the lattices we use is slightly asymmetric, 
30 X 32^. This lattice is close enough to symmetric, however, that the contribution 
of the disconnected term in C2^'{t) turn out to be significantly smaller than our 
statistical errors. 

For any pair of sizes s and s', at sufficiently large values of t and the lattice 
time direction periodicity T, the glueball propagators approach the asymptotic form 

Cf{t) Zf{exp{-m,t) + exp[-mi{T - t)]} (2.9) 

where mo is the energy eigenvalue of the lightest O"*"^ zero-momentum state |0"^"*" >, m2 
is the common energy eigenvalue of the lightest 2++ zero-momentum states \2^^^i > 
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with spin index i = 1,2, and Zq^' and are related to the vacuum state |fi > 
matrix elements 

z^'' = <n\A'{t)\o++ ><o++\A''{t)\n>, 

Zf = ^ < r]|E;(t)|2++,z >< 2++,z|Ef > . (2.10) 

i 

It is not hard to show that, for any value of t, the propagators C-^' (t) are sums 
of terms of the form in Eqs. ( p. 91) with each Z^^' matrix hermitian and non-negative. 
Thus each propagator with the same size smearing at both ends, C-'^lt), is a convex 
function. To extract a glueball mass from each propagator C^^'{t), we search for a 
range of s, s' and t for which the asymptotic form is approached. We then fit the 
asymptotic form to this data set, minimizing the fit's determined from the full 
correlation matrix of the target data set. 

Statistical uncertainties are found for the resulting masses by the bootstrap 
method [^]. From each ensemble of N gauge field configurations, we generate a col- 
lection of 100 bootstrap ensembles. Each bootstrap ensemble consists of N configura- 
tions selected randomly from the underlying N member ensemble, allowing repeats. 
For each of the bootstrap ensembles, masses are determined by repeating our fitting 
method from the beginning but using the fitting interval already chosen for the full 
ensemble. The statistical uncertainty of a mass prediction is taken to be half the 
difference between a mass which is higher than all but 15.9% of the bootstrap masses 
and a mass which is lower than all but 15.9% of the bootstrap masses. In the limit 
of large N, for which the collection of bootstrap masses will approach a gaussian dis- 
tribution, the definition we use for statistical uncertainty approaches the dispersion, 
d, given by (P =< m? > — < m >^. 

Our method of determining masses and error bars is nearly the same as used 
in Ref . ||^ . The sole difference is that in the calculation of Ref . |^ the fitted range of 
s, s' and t is chosen independently on each different bootstrap data set to minimize 
the fit's value of P^r degree of freedom. Although allowing the fitted range of 
s, s' and t to be chosen independently for each bootstrap ensemble appears to be 
logically the best procedure, we were forced to use a fixed fitting range in the present 
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calculation. The statistical fluctuations in glueball propagators for some bootstrap 
data sets were sufficiently large that an automated choice of fitting range would 
be difficult to implement, while choosing fitting ranges by hand for each bootstrap 
ensemble would be too time consuming. 

It is perhaps worth mentioning that the values of glueball masses obtained from 
Eqs. ( p.9|) are independent of the choice of the convergence parameter r. Coulomb 
gauge leaves the time direction link matrices unconstrainted. The integration over 
time direction link matrices in the path integral generates, in the corresponding time 
evolution transfer matrix of the Hilbert space formulation of lattice QCD a 



projection operator onto the gauge invariant sector of Hilbert space. For values of 
r which are large, the states created by our glueball operators will include large 
gauge variant components. These are canceled off by the time evolution projection 
operators to the gauge invariant sector of Hilbert space. The resulting propagators 
of Eqs. ( p.8| ) exhibit only the masses of physical, gauge invariant states. The price 
of this cancellation is an increase in the statistical errors in masses. In practice we 
choose a small enough value of r so that this increase in statistical errors is negligible. 

In comparison to the gauge invarient smearing methods of Refs. H, the 
numerical cost of fixing to Coulomb gauge is rather large. As a results our method 
of constructing glueball propagators is less efficient that those of Refs. |^, An 
advantage of our procedure is that it does a more complete job of eliminating high 
frequency components of the gauge field than do the methods of Refs. |[| and, 
correspondingly, is free of the potential problems discussed in Ref. |T^. A comparison 
of our data with that of Ref. also suggests that, for a fixed ensemble size, our 
method may give somewhat smaller statistical errors in glueball propagators at large 
time separations. 

Our procedure for determining glueball masses from glueball propagators ap- 
pears to be marginally more reliable than two competing methods. We found that 
fitting the full propagator matrix C^^' (t) gives masses which are somewhat more sta- 
ble to variations in ensemble size than either fitting the correlation function Cf*(t) for 
a single optimal smearing size s, or fitting the correlation function J2ss' C^C^' Ci^' (t) 
for a single optimal vector 01- It ^-Iso appears to us that a single fit to the full 
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set of data provides the least biased way of resolving the small disagreements which 
occur among the masses found with various different combinations of smearing size 
for the propagator's source and sink. 

3 MASSES 

Gauge configurations were generated according to the standard Wilson one 
plaquette action using the Cabbibo-Marinari-Okawa algorithm. Table |l] lists the lat- 
tice sizes, (3 = Q/gl values, sweeps skipped between gauge configurations, and number 
of configurations used in the ensembles from which glueball propagators were calcu- 
lated. Each configuration was equilibrated for at least 30000 sweeps before we began 
collecting data. The number of sweeps between configurations on which data was 
collected, however, was not large enough for us to be certain whether or not succes- 
sive data samples were correlated. We therefore checked our bootstrap evaluation of 
error bars to be sure that it was not distorted by the presence of correlations. At (3 of 
5.7, 6.17 and 6.4, our collections of gauge configurations were divided into sequential 
groups of 2, 4, 8 and 16 gauge configurations. Glueball propagators were averaged 
over each bin, and a bootstrap error evaluation was done on masses determined from 
the binned ensembles. Any correlations present in an ensemble of propagators before 
binning become progressively weaker in the ensembles with larger bins. Our error 
estimates showed no significant variation as bins were made larger. A direct calcu- 
lation of the correlation between successive glueball propagators in the ensembles of 
Table ^ also produced no statistically significant evidence for correlations. 

For all the lattices listed in Table the gauge fixing convergence parameter 
r was set to 10~^. The average number of sweeps required for convergence to this 
accuracy ranged from 1525 on the lattice 24^ x 36 at /5 of 5.93, to 2270 on the lattice 
16^ X 32 at /5 of 5.7. 

To determine the range of smearing sizes and time separations from which 
to extract masses for each lattice and spin, we evaluated effective masses m|(t) by 
fitting propagators Cf^it) to Eqs. ( |2.9| ) at separations t and t + 1. In all cases, the 
time direction period T was large enough and the time separation at which our mass 
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signals disappeared into noise small enough that the values we found for m^{t) were 
equivalent to using the simpler definition m^(t)a = hi[C-^(t)/C-^(t + 1)], where a 
is the lattice spacing. At large t each mf{t) approaches the corresponding rrii of 
Eqs. ( p.9|) . It follows from the convexity of C^''^{t) that this approach will be from 
above. A plateau at large t in a graph of Tnflt) identifies the region over which it 
is reasonable to try fitting propagator data to Eqs. (p.9| ). Figures (|l]) - (|) show 
effective masses for the 0"*"^ glueball on all lattices of Table |l| and for the 2"^"^ glueball 
on the three lattices at largest (3. At time separations beyond those shown, we were 
unable to obtain statistically significant effective mass values. At the two lowest f3 
we were unable to obtain enough statistically significant effective mass values for the 
2++ glueball to find a plateau in the mass. In each case, we show the smearing size 
which gives the best signal. All of the data for the 0^"*" glueball shows fairly clear 
effective mass plateaus extending over at least 4 time values. The data for the 2++ 
also shows effective mass plateaus, but not as clearly as the 0^^ data. 

Having found effective mass plateaus for each smearing size, we fit the largest 
set of time slices which could be considered to be in the plateau to Eqs. ( |2.9| ) and 
determined the P^^^ degree of freedom of the best fit. We then removed a time 
slice at small t and repeated the fit. In most cases, either the value of per degree 
of freedom of the best fit decreased or the fitted mass decreased by a statistically 
significant amount, showing that the narrower fitting range provided a better estimate 
of the lowest mass. We repeated the process as long as either the effective mass 
continued to decrease significantly or per degree of freedom decreased. The process 
was stopped when we found either a large increase in per degree of freedom or a 
large increase in the fitted mass's statistical uncertainty. We adopted the time window 
before either a large increase in or a large increase in statistical uncertainty as the 
best fitting range. In about half of the cases we considered, a fit to the window with 
one fewer time slice than the best choice showed a larger value of mass than the best 
choice. The convexity of the propagators with the same smearing size at both ends 
implies this increase can only occur as a result of statistical noise. The occurrence 
of an increase in fitted mass about half of the time suggests that, within statistical 
errors, we obtained the correct lowest mass in each case. Tables (0) and (|^) show the 
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masses and values of per degree of freedom obtained in the process of choosing 
fitting ranges for the data in Figures (P - (||). For each sequence of fits in Tables (0) 
and (H), the fitting range next to last was selected as the best choice. The vertical 
lines in Figures ([^) - show the final fitting range and the horizontal lines show 
the mass found for this range. In some cases the fitting range extends beyond the 
last time slice for which an effective mass is given since, by convention, we assign 
each effective mass to the smaller of the two time values from which it has been 
determined. Tables @, @, and @ show the masses found by fitting, on the best 
fitting ranges, propagators with the same size smearing at both ends. For each (3 and 
spin, the optimal fitting range determined according to our criteria turned out to the 
be same for each smearing size for which we had a reliable signal. This was not put 
in as a constraint, is slightly surprising, and is probably a statistical accident. 

Final mass values for each particle were determined by making a combined 
fit, as described earlier, to data for the full set of different smearing sizes for which 
we had reliable masses. The results of the final fits are shown in Tables and ®. 
For the 2++ glueball, at P of 6.17, the combined data set dictated by our selection 
rule consisted of smearing sizes 4, 5 and 6 over the range from time separation 2 to 
time separation 5. This fit gave an unacceptable per degree of freedom of 10.6. 
Removing the data at time separation 2 and fitting only from separations 3 to 5 
reduces of 1-2 and gives the mass shown in Table (|^). 

In addition to our calculations of propagators and masses for the 0"'""'" and 2++ 
glueballs, we also calculated propagators for the 1~^~ glueball. The effective mass 
plateaus in these propagators were weak, however, and the masses which could be 
extracted were close to the expected masses for torelon pairs. Thus the reliability of 
the predicted 1"*"" masses was not clear, and we do not report these results here. 

4 INFINITE VOLUME, CONTINUUM LIMITS 

The quantity z = mo++L, where L is the lattice period, for the masses in Table 
(1^ ranges from 12.45, using the smaller period of 30 on the lattice 30 x 32^ x 40 at 
p of 6.4, to 19.46, on the lattice 24^ x 36 at /5 of 5.93. Masses calculated in a large 
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finite volume approach their values in infinite volume according to [|T2 



m,{z) = m,(oo)[l - ^,f^^U!fl]. (4.1) 



For mo the coefficient has been estimated |13| to be 190 ± 70. Thus our values 



for mo should differ from their infinite volume values by less than 0.05%. For m2 it 



appears from data in Ref. |T3| that g2 < 2000. The range of z we use then gives m2 
which agree with infinite volume values to better than 0.4%. Thus all of the masses 
shown in Tables (|^) and differ from their values in infinite volume by amounts 
which are negligible in comparision to the statistical uncertainties. 

In preparation for evaluating the continuum limit of uiq and m2 in physical 
units, the values we calculated for m^a in lattice units were rescaled by isl^a in 

' ' MS 

lattice units. According to mean-field improved lattice perturbation theory [jl4|, |15 

1 < TrU/3 > 



9hs 9' 



+ 0.025, (4.2) 



where g = Q/(3 is the lattice coupling constant, and < TrU > is the average value of 
a plaquette. Then with ajjg = gj^/Aii, the solution to the two loop Callan-Symanzik 
equation for aj^{iT/a) with zero fiavors of quark vacuum polarization gives 



bo = ^, (4.3) 



51 

22^' 



Values of aj^g and A^^a determined from Eqs. ( |4.2| ) and ( ^^ ) are given in Table 
(^ . Values of m^ / A^^ determined from Tables (|^) , (P) and (^) are shown in Figures 
(H) and (|10D- In Ref. 0] it was found that the rho mass ratio mp/A^^ is nearly 
independent of (3 and within 2% of its continuum limit once (3 > 5.7. Thus graphs 
of our glueball masses scaled by the rho mass, mi/rrip, would differ from Figures (^ 
and (ll^) only by a constant rescaling of the vertical axis. 
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Since the leading irrelevant term in the Wilson action for pure gauge fields is 
O(a^), we expect that for sufficiently small lattice spacing mass ratios will approach 
their continuum limits according to 



with a- independent coefficients q. Fits of ttiq to Eq. ([4 .41) using five data points, using 
the four data points with the smallest a and using the three data points with the 



progressively decreases as data points with larger a are eliminated, it appears that the 
fit using three data points gives the most reliable extrapolation to the continuum limit. 
It is interesting to notice that a quadratic fit to all five data points gives essentially 
the same continuum limit and same P^r degree of freedom as the linear fit at the 
three smallest points. For 1712 we have data at only the three smallest values of a. A 
fit of Eq. ( [4.4|) to these, giving a fairly small is also shown in Table ([T0|). The three 
point fits for mo and m2 are shown in Figures @ and (0), respectively. The vertical 
lines at a = in each figure are the statistical uncertainties in the extrapolated mass 
values. 

To obtain masses in physical units from the continuum limit ratios mj/A^^, 
we use the continum value of 245.0 ± 9.2 MeV for found in Ref. from the 
continuum limit of mp/A^^ combined with the observed value of rup. For mo we get 
1340 ± 160 MeV, and for m2 we find 1900 ± 320 MeV, with statistical uncertainties 
combining the uncertainties in mJK^^ and in A-^. 

° MS MS 

An alternate way to find the continuum limit of lattice values of mjO is to divide 
each by the square root of the string tension in lattice units, \fKa, and extrapolate 
mi/\fK to its continuum limit. The continuum limit of \fK, however, is not a directly 
observable quantity. Its continuum value can be determined by extrapolating \fK /irtp 
to its continuum limit, then using the observed value of rrtp. The resulting continuum 
prediction for \fK inherits statistical errors both from lattice calculations of \fK and 
from lattice calculations of nip. The continuum value which we have used for A^^, 




MS 



MS 



(4.4) 



smallest a are listed in Table (p!0|). Since Eq. ( [4.4| ) is an asymptotic form approached 
at small lattice spacing a, and since the per degree of freedom of fits to Eq. ( [4.4| ) 
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on the other hand, has statistical errors only from lattice calculations of nip. Thus 
the method we have adopted should give a smaller uncertainty for continuum mass 
predictions in physical units. The central value of the answers obtained extrapolating 
mi/\fK do not appear to be significantly different from our present answers. For 



example, interpolating pubhshed values of yKa to our values of /5, using these 
to find the continuum limit of mQ/\fK, and assuming the popular value of 440 MeV 
for continuum V^, gives a continuum mo of about 1310 MeV. 

As mentioned in Sect. (0), continuum limit predictions for mo of 1550±50 MeV 
and for m2 of 2270 ± 100 MeV were reported recently in Ref. These numbers were 
obtained by taking mo/V^ and 7712/ \/K at (3 of 6.4 as continuum values. If the 
fit in Ref. of vtiqI^/K to a linear function of Ka^ is actually extrapolated to the 
continuum, the central value of mo becomes 1600 MeV. This number differs from our 
result by about 1.6 standard deviations. It appears to us, however, that the central 
values of both predictions quoted in Ref. may be somewhat higher and the error 
bars on both predictions somewhat smaller than the data of Ref. |^ itself suggests. 

For mo and m2 at /5 of 6.4, Ref. ||^ takes the effective masses found from the 
0"*"+ and 2++ propagators between time separations 2 and 3. From the effective mass 
tables given in Ref. |^ it seems unclear how close the effective masses at these time 
separations are to their asymptotic large time separation values and how reliable the 
error estimates on these effective masses may be. The central values of effective masses 
found from the 0"*"+ and 2++ propagators fall monotonically at time separations larger 
than 2. Both effective masses between time separations of 4 and 5, for example, are 
about 10% lower than corresponding numbers between separations 2 and 3 and have 
error bars which are larger than the corresponding errors by factors of more than 2. It 
also seems to us uncertain whether the mo value which Ref. takes from Ref. for 
/5 of 6.2 is a reliable estimate of the asymptotic effective mass at large time separation. 
Our calculation of mo at (i of 6.17, which should lie above mo at 6.2, is actually lower 
by about 10%. The data of Ref. for mo at of 5.9 is consistent with our fit of 
mo to a linear function of [A^^a]^, while m2 is again above the extrapolation of our 
data. 

The results of Refs. [0 and become closer to ours if mo and m2 are taken 
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from effective masses at larger time separations approximating the ranges over which 
our mass fits were done. Combining effective masses between time separations 2 
and 3 at /3 of 5.9 and 6.0 0, time separations 3 and 4 at /3 of 6.2 and time 
separations 4 and 5 at /5 of 6.4 p, we fit mj/A^^ to hnear functions of [A^^a]^. The 
continuum mo/A|^ becomes 6.11 ± 0.52 and m2/A^^ becomes 7.75 ± 0.89. With 
^jjg of 245.0 ± 9.2 MeV, the continuum prediction for uiq is 1500 ± 140 MeV and for 
7712 is 1900 ± 230 MeV. These results are both less than one standard deviation from 
our corresponding predictions. A fit to this set of data combined with ours in the 
region (3 > 5.9 gives a continuum limit for mo/A^^ of 5.89 ± 0.38 and for 7712/ 
of 7.77 ± 0.73. These values yield mo of 1440 ± 110 MeV and 7712 of 1900 ± 190 MeV. 

We would like to thank Frank Butler for writing some of the analysis software 
which we used, and Mike Cassera, Molly Elliott, Dave George, Chi Chai Huang and 
Ed Nowicki for their work on GFll. 
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lattice 


/3 


skip 


count 


16=^ X 24 


5.70 


400 


4039 


20^ X 30 


5.83 


400 


4002 


243 X 36 


5.93 


400 


4004 


30 X 32^ X 40 


6.17 


400 


2005 


30 X 322 X 40 


6.40 


400 


2002 



Table 1: Configurations analyzed. 





size 


range 


mass 




5.70 


2 


1-5 


0.958 ±0.023 


0.56 


5.70 


2 


2-5 


0.896 ± 0.058 


0.20 


5.70 


2 


3-5 


0.978 ±0.176 


0.03 


5.83 


2 


1-5 


0.896 ±0.017 


0.93 


5.83 


2 


2-5 


0.845 ± 0.039 


0.62 


5.83 


2 


3-5 


0.798 ± 0.076 


0.92 


5.93 


2 


1-5 


0.856 ±0.015 


1.15 


5.93 


2 


2-5 


0.833 ±0.033 


1.54 


5.93 


2 


3-5 


0.871 ±0.085 


2.78 


6.17 


4 


2-7 


0.564 ±0.028 


1.55 


6.17 


4 


3-7 


0.482 ± 0.035 


0.20 


6.17 


4 


4-7 


0.461 ±0.058 


0.20 


6.40 


5 


3-9 


0.444 ± 0.029 


1.12 


6.40 


5 


4-9 


0.387 ±0.035 


0.34 


6.40 


5 


5-9 


0.399 ±0.050 


0.43 



Table 2: Masses, in units of 1/a, obtained for O"*""*" glueballs over various fitting ranges. 
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p 


size 


range 


mass 


2 

X 


5.93 


4 


1-4 


1.22S ± 0.050 


0.03 


o.yo 


4 


O /I 


1 1 n7 _L n 1 QQ 
i.iy/ ± U.ioo 


n no 


5.93 


4 


3-4 


1.288 ±0.766 




6.17 


5 


1-5 


0.871 ± 0.029 


0.22 


6.17 


5 


2-5 


0.852 ± 0.061 


0.28 


6.17 


5 


3-5 


0.835 ± 0.137 


0.55 


6.40 


5 


2-8 


0.579 ± 0.023 


1.25 


6.40 


5 


3-8 


0.632 ± 0.053 


1.09 


6.40 


5 


4-8 


0.522 ± 0.059 


0.49 


6.40 


5 


5-8 


0.641 ± 0.128 


0.04 



Table 3: Masses, in units oil /a, obtained for 2++ glueballs over various fitting ranges. 



p 


size 


range 


mass 


x' 


5.70 





2-5 


0.985 ±0.141 


0.01 


5.70 


1 


2-5 


0.948 ± 0.064 


0.54 


5.70 


2 


2-5 


0.896 ±0.058 


0.20 


5.70 


3 


2-5 


0.842 ± 0.072 


0.57 


5.70 


4 


2-5 


0.797 ±0.103 


0.87 


5.83 


1 


2-5 


0.889 ±0.056 


2.02 


5.83 


2 


2-5 


0.845 ± 0.039 


0.62 


5.83 


3 


2-5 


0.840 ± 0.047 


1.30 


5.83 


4 


2-5 


0.845 ± 0.075 


1.36 


5.83 


5 


2-5 


0.834 ±0.116 


0.79 


5.83 


6 


2-5 


0.798 ±0.105 


1.00 



Table 4: 0+"*" glueball masses, in units of 1/a, fitted to optimal ranges for each size 
of smearing. 
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/9 


size 


range 


mass 




5.93 


1 


2-5 


0.905 ±0.053 


0.13 


5.93 


2 


2-5 


0.833 ±0.031 


1.54 


5.93 


3 


2-5 


0.796 ±0.030 


3.19 


5.93 


4 


2-5 


0.781 ±0.043 


3.98 


5.93 


5 


2-5 


0.746 ± 0.060 


2.41 


5.93 


6 


2-5 


0.724 ± 0.090 


1.52 


6.17 


2 


3-7 


0.504 ± 0.055 


0.83 


6.17 


3 


3-7 


0.505 ± 0.035 


0.37 


6.17 


4 


3-7 


0.482 ± 0.035 


0.20 


6.17 


5 


3-7 


0.473 ± 0.041 


0.04 


6.17 


6 


3-7 


0.495 ± 0.060 


0.30 


6.40 


3 


4-9 


0.392 ± 0.043 


0.38 


6.40 


4 


4-9 


0.389 ± 0.037 


0.27 


6.40 


5 


4-9 


0.387 ±0.035 


0.34 


6.40 


6 


4-9 


0.392 ± 0.034 


0.52 


6.40 


7 


4-9 


0.400 ± 0.036 


0.66 


6.40 


8 


4-9 


0.408 ± 0.042 


0.77 



Table 5: 0++ glueball masses, in units of 1/a, fitted to optimal ranges for each size 
of smearing. 
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size 


range 


mass 




5.93 


3 


2-4 


1.131 ±0.109 


0.11 


5.93 


4 


2-4 


1.197 ±0.188 


0.03 


5.93 


5 


2-4 


1.464 ±0.405 


0.19 


6.17 


4 


2-5 


0.842 ±0.059 


0.18 


6.17 


5 


2-5 


0.852 ±0.062 


0.28 


6.17 


6 


2-5 


0.887 ±0.080 


0.22 


6.40 


4 


4-8 


0.514 ±0.066 


0.41 


6.40 


5 


4-8 


0.522 ±0.059 


0.49 


6.40 


6 


4-8 


0.510 ±0.062 


0.81 



Table 6: 2"'""'" glueball masses, in units of 1/a, fitted to optimal ranges for each size 
of smearing. 





sizes 


range 


mass 




5.70 


0-4 


2-5 


0.928 ±0.058 


1.35 


5.83 


1-6 


2-5 


0.858 ± 0.043 


1.80 


5.93 


1-6 


2-5 


0.811 ±0.033 


2.70 


6.17 


2-6 


3-7 


0.489 ±0.031 


1.29 


6.40 


3-8 


4-9 


0.415 ± 0.043 


1.68 



Table 7: 0++ glueball masses, in units of 1/a, obtained from fits to propagator ma- 
trices for sets of several smearing sizes. 
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/3 


sizes 


range 


mass 




5.93 


3-5 


2-4 


1.144 ±0.107 


2.58 


6.17 


4-6 


3-5 


0.816 ±0.119 


1.17 


6.40 


4-6 


4-8 


0.504 ±0.061 


2.22 



Table 8: 2++ glueball masses, in units of 1/a, obtained from fits to propagator ma- 
trices for sets of several smearing sizes. 



lattice 


/5 






16^ X 24 


5.70 


0.1456 


0.1661 


20^ X 30 


5.83 


0.1370 


0.1329 


24=^ X 36 


5.93 


0.1318 


0.1144 


30 X 32^ X 40 


6.17 


0.1218 


0.08265 


30 X 322 X 40 


6.40 


0.1141 


0.06177 



Table 9: Values of a-jjg and A^g^^a. Statistical uncertainties are all smaller than 1 in 
the last decimal place. 



jPC 


fit 


points 






0++ 


linear 


3 


5.46 ± 0.64 


2.36 


0++ 


linear 


4 


6.19 ±0.50 


2.89 


0++ 


linear 


5 


6.97 ±0.29 


3.43 


0++ 


quadratic 


5 


5.37 ±0.81 


2.26 


2++ 


linear 


3 


7.76 ± 1.26 


0.47 



Table 10: Extrapolations of nii/ Aj^ to the continuum limit using either linear or 
quadratic functions of [Aj^^a]^. 
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Figure 1: Effective masses and fitted mass for the 0+"^ glueball on a 16^ x 24 lattice 
at /3 = 5.7 for smearing size 2. 
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Figure 2: Effective masses and fitted mass for tfie 0''"'*' glueball on a 20^ x 30 lattice 
a,t P — 5.83 for smearing size 2. 
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Figure 3: Effective masses and fitted mass for the 0+"^ glueball on a 24^ x 36 lattice 
a,t P — 5.93 for smearing size 2. 
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Figure 4: Effective masses and fitted mass for tfie O'*"'" glueball witfi smearing size 4 
on a 30 X 32^ x 40 lattice at /3 = 6.17. 
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Figure 5: Effective masses and fitted mass for tfie O'*"'" glueball witfi smearing size 5 
on a 30 X 32^ x 40 lattice at /3 = 6.40. 
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Figure 6: Effective masses and fitted mass for the 2++ glueball with smearing size 4 
on a 243 X 36 lattice at /3 = 5.93. 
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Figure 7: Effective masses and fitted mass for the 2++ glueball with smearing size 5 
on a 30 X 32^ x 40 lattice at /3 = 6.17. 
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Figure 8: Effective masses and fitted mass for tfie 2+"^ glueball witfi smearing size 5 
on a 30 X 32^ x 40 lattice at /3 = 6.40. 
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